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ABSTRACT 



The present invention comprises a system and method for a 
GSM receiver to perform channel estimation under the 
assumption that the analog-to-digital (A/D) converter is 
free-running. A search process is employed, whereby the 
known GSM training signal is sampled with a plurality of 
phases which vary incrementally from zero to a full symbol 
period. For each phase a sample set (of the training signal) 
is generated. Furthermore, each sample set is used, together 
with the received samples, in a cross-correlation procedure 
to obtain a candidate impulse response vector for the trans- 
mission channel. Thus a plurality of candidate impulse 
response vectors are produced, one of which will be selected 
to be the optimal representative for the transmission chan- 
nel. The optimal impulse response vector is selected as 
follows. Each impulse response vector is convolved with the 
corresponding sample set, thus producing an estimated 
output vector. The estimated output vectors are compared 
with the vector of received samples (from the free-running 
A/D converter), and the estimated output vector which is 
closest in the Euclidean (square-norm) sense designates the 
impulse response vector which is optimal. 

53 Claims, 7 Drawing Sheets 
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SYSTEM AND METHOD FOR ESTIMATING 
A SET OF PARAMETERS FOR A 
TRANSMISSION CHANNEL IN A 
COMMUNICATION SYSTEM 

5 

CROSS REFERENCE TO RELATED 
APPLlCAnON(S) 

This application claims benefit of priority of provisional 
application Serial No.: 60/040,195 filed Mar. 4, 1997, titled 
Robust Channel Estimator for GSM Applications. 

FIELD OF THE INVENTION 

This invention is related to the field of channel estimation 
in the context of cellular radio communication, and particu- 15 
larly to a channel estimation algorithm for a receiver struc- 
ture in the GSM cellular radio system. 

DESCRIPTION OF THE RELATED ART 

20 

In the GSM cellular radio system, base stations and 
mobile units transmit and receive digitized voice data in the 
900 MHz band. Tht 890-915 MHz band is reserved for 
mobile -to-base transmissions (up link), and the 935-960 
MHz band is reserved for base-to-mobile transmissions 
(down link). These bands are divided into channels which 
are 200 kHz wide. A pair of channels, an up link channel and 
a down link channel, separated by 45 MHz is described by 
an Absolute Radio Frequency Channel Number (ARFCN), 
In normal communications each mobile unit is assigned an 
ARFCN, asserting its transmissions over the up link 
channel, and receiving base station transmissions over the 
down link channel. 

Each base station is capable of operating over several 
ARFCNs, Furthermore, on each ARFCN the base station 35 
may be in contact with up to eight mobile units. The mobile 
units share an ARFCN by a scheme of Time-Division 
Multiple Access (TDMA), wherein communications over 
the up link and down link channels are divided into frames, 
and each frame is further sub-divided into eight time -slots. 40 
Each of the eight mobile units is assigned a time-slot number 
and is expected to confine its transmissions on the up link 
channel to the corresponding time-slot, and to receive down- 
fink transmissions in the corresponding time-slot. Thus the 
GSM transmitter (base or mobile) must assemble speech 45 
data into bursts (of length 148 bits). When asserting the burst 
during a time-slot, the transmitter shapes the power enve- 
lope of the transmitter RF signal in conformity with a GSM 
power-ramping time-mask [ETSI GSM 05.05, "European 
Digital Cellular Telecommunication System (Phase 2), 50 
Radio Transmission and Reception", p.34, September 1995]. 
The power envelope is used by the receiver (base or mobile) 
to detect the start of the burst and to synchronize its A/D 
sampling. 

The GSM transmitter uses the Gaussian Minimum-Shift 55 
Keying (GMSK) modulation technique, in which the binary 
speech data is embedded onto a baseband signal that modu- 
lates the carrier with a bit rate of 270833 bits per second. The 
modulated carrier is transmitted into the atmosphere as an 
RF signal. The receiver accepts the RF signal from the 60 
atmosphere, and down-converts this RF signal to produce a 
baseband signal. Due to the effects of the atmosphere or 
transmission channel, the received baseband signal is a 
distorted version of the transmitted baseband signal. The 
distortions include atmospheric channel fading, multipath, 65 
noise, and the distortions imposed by the RF stages of the 
transmitter and receiver. 
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A typical GSM receiver structure includes: an RF stage 
which converts the received RF signal into a baseband 
signal; an analog-to-digital (A/D) converter which samples 
the baseband signal at a rate of 270833 samples per second, 
i.e. once per bit interval, to produce a sequence of received 
samples; a channel estimator which estimates parameters 
describing the transmission channel on a burst by burst 
basis; and a Viterbi equalizer which uses the channel param- 
eters to demodulate the binary speech data from the A/D 
samples. 

In the 900 MHz band transmissions are subject to channel 
fading and multipath, which results in significant inter- 
symbol interference (ISI). Furthermore, the RF stages of the 
transmitter and the receiver are non-ideal, and exhibit their 
own transfer characteristics. The channel estimator is 
designed to estimate parameters which describe the trans- 
mission charmel, i.e. estimate the transfer function of the 
transmission charmel. In this appUcafion, the transmission 
channel is interpreted to include one or more (preferably aU) 
of the signal stages between the transmitted baseband and 
received baseband signals. Thus the transmission channel 
includes one or more of the following: the RF stage of the 
transmitter, the atmospheric medium, and the RF stage of the 
receiver. The channel estimator assumes that the transmis- 
sion channel can be modeled as a discrete-time finear 
system, and calculates the impulse response of this linear 
model. 

To facilitate the estimation procedure for the transmission 
channel impulse response, each transmitter burst contains 
known training data. In particular, every burst, which com- 
prises 148 successive bits, contains a training sequence of 
length 26 bits. There are eight different training sequences 
used in GSM communications, all of which were chosen 
because of their nearly optimal auto-correlation behavior. 
The channel estimator uses the samples of the received burst 
and the known training sequence in its estimation algorithm. 
The Viterbi equalizer uses the transmission channel impulse 
response or transfer function calculated by the channel 
estimator to process the received samples and compensate 
for the ISI. In doing so, the \^terbi equalizer demodulates 
the encoded binary speech data. 

In order to fit a discrete linear model to the transmission 
channel, the channel estimator must provide the finear 
model with an input and output sequence. The input 
sequence comprises samples of the transmitted baseband 
signal, and the output sequence is quite naturaUy determined 
by the received digital samples output from the A/D con- 
verter. In prior art systems it has been the practice to control 
the A/D converter so that samples of the received baseband 
signal are obtained at favorable locations within each bit 
interval; i.e. it has been necessary to control the phase of the 
A/D converter (not the sampling rale). The favorable loca- 
tions are determined ultimately by the bandwidth of the 
Gaussian pulse-shaping filter used in the GMSK modula- 
tion. By obtaining a sample 2.5 symbol periods after the 
assertion of each binary digit, the ISI is minimized [see page 
172 of: Wireless Digital Communications: Modulation and 
Spread Spectrum Applications by Dr. Kamilo Feher, ISBN 
0-13-098617-8]. Such practice has required prior art GSM 
receivers to include circuitry for recovering bit synchroni- 
zation. Furthermore the A/D converters in prior art receivers 
have required a control input and extra logic to allow for 
sample liming adjustment. 

Thus a system and method is needed whereby the channel 
estimation can be performed successfully without control- 
fing A/D sample timing. Such a method would allow the 
receiver structure to be simplified; and the A/D converter, 
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not neediag the extra control input and logic, would con- corresponding sample set, tfaus producing an estimated 

sume less power and occupy less semiconductor die area. output vector. The estimated output vectors are compared 

with the vector of received samples, and the estimated 

SUMMARY OF THE IiWENTION output vector which is closest, preferably in the EucUdean 
Hie present invention comprises a system and method for ^ (square-nomi) sense, designates the impulse response vector 

performing channel estimation in a GSM receiver. The "^^^^^ ^ optmial. This impulse response vector, which best 

channel estimation operates under the assumption that the describes the transmission channel, can then be used to 

analog-to-digital (A/D) converter is free-running. Thus, demodulate the binary data embedded m the received 

samples of the received baseband signal are obtained with samples. 

unknown (but constant) phase. The received baseband signal nP^rRTimnM HP thp hr awtmh^ 

is a distorted version of a transmitted baseband signal, due ^^^^^ DESCRIPTION OF THE DRAWINGS 

to the effects of the transmission channel. During the train- a better understanding of the present invention can be 

ing phase of every GSM signal burst, the transmitted base- obtained when the following detailed description of the 

band signal is completely known to the receiver: this known preferred embodiment is considered in conjunction with the 

signal is called the training signal. Since the phase used to following drawings, in which: 

obtain the received samples is not controlled, it is not loiown ^ ^ ^^^'^^ ^-^^^^^ ^gM transmission 

aprion which samples of the trammg signal would best system 

correspond with the received samples. In other words, since • j- . r ^^.^m ^ 

*u u a- * f .L -J 1 • .1 • FIG. 2 is a diagram of the structure of a GSM data burst. 

the phase onset of the received samples is not known, it is ^ 

not known which samples of the transmitted training signal ^ FIG. 3 is a block diagram of a GMSK Modulator, 

should be used in correlation with the received samples. FIG. 4 is an illustration of the means by which sampling 

In the system and method of the present invention a search ^ performed on the known training signal, to generate a 

process is employed. First, the training signal is sampled plurality of sample sets. 

with a plurality of phases which vary incrementally from ^ FIG. 5 (referring collectively to FIGS. 5 A, 5B, and 5C) 
zero to a full symbol period. For each phase a sample set of presents flowcharts which illustrate the method of the 
the training signal is generated. Furthermore, each sample present invention. FIG. 5A presents the invention in over- 
set is used, together with the received samples, in a cross- view. FIG. SB presents the preferred embodiment of the 
correlation procedure to obtain a candidate impulse response channel estimation process. FIG. 5C presents an alternate 
vector for the transmission channel. Thus, a plurality of embodiment of the channel estimation process, 
candidate impulse response vectors are produced, one for pic. 6 is a flowchart which is used to motivate a discus- 
each of the phase offsets. Finally, one of the impulse sion of possible alternate embodiments of the parameter 
response vectors is selected to be the optimal representative e.stimation process, 

for the transmission channel. j -3 ^ ^^^^^ ^-^^^^^ ^ ^-^^y ^^^^^ processing 
The cross-correlation procedure operates on a sample set 35 system according to the preferred embodiment of the present 

and the received samples to produce a candidate impulse invention. 

response vector; one candidate itnpulse response vector is ^^i^ j^^^^jj^^ ^ susceptible to various modifica- 

produced for each sample set. In the preferred embodiment, ^^^^ alternative forms specific embodiments are shown 

the cross-correlauon procedure mvolves performing two ^ ^^^^ ^^^^^ ^ ^^^j^ ^^ 

cross-correlauons: a cross-correlaUon between the received a^scvib^d in detail. It should be understood however, that 

samples and the known bmary data (embedded in the ^^^^ j^.^jj^ ^ description thereto are not intended to 

traming spal), and a cross-correlation between the sample ^^-^ ^^^^^^^^ ,^ ^^^^ particular form disclosed. But on 

set and the same binary data. These two «,rrelations are ^ ^ -^^^^^^^ ^ ^^^^^ ^, ^^^^^^^^^^^ 

connected by an underlying Imear system. By solving this equivalents and alternatives following within the spirit and 

underlying Imear system, a candidate unpulsc response ,5 j invention as defined by the appended 

vector IS produced. In particular, the hnear system is solved claims 
by matrix inversion. A vector is composed from the results 

of the first cross-correlation, and a matrix is composed from DETAILED DESCRIPTION OF THE 

the results of the second cross-correlation. The matrix is PREFERRED EMBODIMENT 
inverted, and the inverted matrix is multiplied by the vector, jq 

The product specifies the candidate impulse response vector. Referring now to FIG. 1, a block diagram of a GSM 

Applying this correlation procedure for each sample set, transmission system is shown. The present invention is 

a plurality of candidate impulse response vectors is gener- preferably used in a GSM system, and the invention is 

ated. The first cross-correlations between the received described below with reference to a GSM system. However, 
samples and the known binary data need be performed only 55 ^ ^^^^ ^^^^ present invention can be used iri other 

once for a received signal burst. Furthermore since the communication systems and/or with other communication 

second cross-correlations involve the known sample sets and standards, as desired. 

the known binary data, the inverted correlation matrices may The GSM transmission system comprises a GSM trans- 
be calculated and stored pre-operationally (at compile time mitter 100, a Radio Channel 140, and a GSM receiver 180. 
or design time, before operation of the system). Thus the The GSM transmitter 100 comprises an Encoder 110, a Burst 
calculation of the impulse response vectors involves per- Assembler 115, a GMSK Modulator 120, and an RF Trans- 
forming the first cross-correlation (of received sample and mitter 130. The GSM Receiver 180 comprises an RF 
binary data) once per received signal burst, and multiplying Receiver 150, an Analog-to-Digital (A/D) Converter 155, a 
the resulting correlation vector by. a plurality of stored Channel Estimator 160, and a Viterbi Eqtializer 170. 
mverse matnces. 65 fhe Encoder 110 receives a stream of sampled voice data 
The optimal impulse response vector is selected as fol- and performs speech and channel coding as well as encryp- 
lows. Each impulse response vector is convolved with the tion. The resulting bit stream is assembled into bursts of 
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leogth 148 bits by the Burst Assembler 115. The assembled because of their nearly optimal aulo-correlalion behavior, 
bursts are modulated using Gaussian Minimum-Shift Key- The bursts arc modulated in the GMSK Modulator 120 using 
ing in the GMSK Modulator 120. The rcsuUing baseband Gaussian Minimum-Shift Keying (GMSK). 

signal is imaged onto an RF carrier and transmitted by the n r * » rtr- ^ i_i i j- c*u r-nxov 

n^T- •** * J nr- • i • Referrmg now to FIG. 3, a block diagram of the GMSK 

RF Transmitter 130. The transmitted RF signal expcnenccs s^,,, . ^ , , 

ahostofdistortionphenomenonasitpropagateslhroughthe Modulator 120 is shown. The GMSK Modulator 120 
atmosphere, all of which arc summarized by the functional ^^^^ Assembler 115 a sequence of 148 

block denoted Radio Channel 140. The distortions include successive bits denoted {ao,ai,a^a3, . . . ^^47} where each 
channel fading, multipath. and noise. The RF Receiver 150 bit a;, takes the value 0 or 1. Before moduladon however, the 
accepts the RF signal from the ataiosphere, and down- jo symmetrized and differentially encoded in Symme- 

converts this RF signal to generate a baseband signal. The tnze and Differential Encode Unit 310. By applying the 
received baseband signal is a distorted version of the trans- transformation b^=2ajt-l a symmetric bit stream {bojbjjb^, 
mitted baseband signal. bg, . . . ^b^^^} is generated, Le. the symmeUized bits take the 

The A/D Converter 155 samples the received baseband values 1 or -1. The symmetric bit stream is differentially 
signal at a rate of 270833 bits per second, and supplies the 15 encoded according to the equation 
samples to the Channel Estimator 160 and Viterbi Equalizer 
170. According to the present invention, the A/D converter 
155 is not required to receive control signals to control the d^^i^kK-^ 
sample timing. Therefore, the A/D converter 155, not need- 
ing extra control inputs and logic, can be designed so as to 20 with b.j defined to be -1, thereby generating bit stream 
consume less power and occupy less die area than A/D {do,di,d^d3, , . . , d^^-f}. 

converters in prior art GSM receivers. The A/D converter differentially encoded bit stream is then formed into 

155 samples at a constant rate, but at an unknown phase ^ bipolar NRZ signal by pulse-generation filter 320. This is 
offset within each bit interval. accomplished by applying the bit stream {do,di,d2,d3, . , . , 

The Channel Estimator 160 uses the baseband signal 25 dj^^} to a filter whose impulse response p(t) is a unit 
samples to calculate, on a burst by burst basis, the impulse amplitude rectangular pulse with turn-on interval O^t^T, 
response or transfer function of the transmission channel, where T is the bit encoding interval, also called symbol 
which includes one or more of the RF Transmitter 130, the interval, equal to 3.692 /cec. The differential bit stream is 
Radio Channel 140, and the RF Receiver 150. The Channel suppUed to the pulse-generation filter as a series of impulses 
Estimator 160 includes an improved system and method 30 
according to the present invention for estimating the channel 147 
parameters, i.e. the transmission channel impulse response ^rft(5(f-Ar), 
or transfer function, even though the phase offset of the 
received signals is unknown. 

The Channel Estimator 160 provides the impulse response, 35 . . . , w ^n. *u xmr? • 1 

, 11 J *u u 1 *■ . u 1 * * where o(t) is the dirac delta function. Thus the NRZ signal 

also called the channel esUmate or channel parameters, to u * ^ 

the Viterbi Equalizer 170. The Viterbi EquaUzer 170 uses the represented as 

channel estimate in a Maximum Likelihood Sequence Esti- 
mation algorithm to demodulate the data bits from the ^^t^pit ~ kry 
received baseband samples. 40 * 
GSM Transmitter Background 

TTie GSM Transmi^ion system mustrated in FIG. 1 ^ signal is then smoothed with a Gaussian 

applies to transmission both from mobile umt to base station ,^ j ^^^^ ^^^^ j ^ 
and from base station to mobile vinit. Tlie signal processing Gaussian function defined by 
structure is similar in both directions. For a complete expo- 45 

sition of the invention, background on the operation of a 

GSM transmitter is deemed appropriate. ^^^^ = B } 

The functional elements of the GSM Transmission system V 1^2 

will be described in detail. For those functional elements 
which are not relevant to the present invention, the corre- 50 

spending description has been omitted for simpUcity. with time-bandwidth product BT equal to 0.3. Thus function 

The Burst Assembler 115 receives a bit stream from the g(0 has even symmetry about time t=2T. Furthermore, 

Encoder HO, and packages the bit stream into bursts of outside the time-interval 0it^4T, the function g(t) is essen- 

length 148 bits. Since GSM specifies Time Division Mul- tidily zero. The smoothed NRZ signal can be represented as 

tiple Access (TDMA) signaling, i.e. multiple mobile units 55 

share a common frequency channel by a scheme of time- ^i/tKr-AH 

domain multiplexing, GSM transmitters must transmit data t 

by means of time-limited signal bursts. 

Referring now to FIG. 2, a diagram of the burst structure 

is shown. A burst comprises two sub-blocks of encrypted 60 "^^^'^ ^0 ^ the convolution of p(t) and g(t). 

data 200 and 210, each of length 58 bits. The two sub-blocks The smoothed NRZ signal is then applied to an integrator 

arc separated by a sequence of 26 training bits 220. Fur- 340. The output of the integrator can be represented as 

thermore the sub-blocks are bounded by tail bit sequences 

230 and 240, each of length 3 bits. The training sequence is Jfc7> 

included in every burst to facHitate channel estimation by 65 l 

the receiver. In GSM communications there are eight dif- 
ferent training sequences, all of which have been chosen 
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where q(t) is the rcsuli of applying r(i) to an integrator, Le. The baseband signal y(l) is sampled by the A/D converter 

155 with sampling rate 270833 samples per second, to 

_ r/<5)rfj produce a sequence of received samples {yo,yi,y2»y3 • • • » 

Jo Yw?}- ^oie that the A/D converter 155 obtains one sample 

5 per bit interval (also called symbol interval). 

The function q(t) rises smoothly from zero to one in the time The Channel Estimator 160 according to the Present Invcn- 

interval from zero to 5T, and has odd symmetry about time tion 

t=2.5T. Operation of the Channel Estimator 160 of the present 

The output of the integrator is multiphed by n/2 to ^^^^^^^ ^ described in detail. According to the 

generate phase signal p^^^^^ invention, the Channel Estimator 160 uses the 

received samples and its apriori knowledge of the training 

^0= ~ydi,gu~kT^. sequence bits to estimate the impulse response (or transfer 

^t=o function) of the transmission channel. In the present 

15 disclosure, the term "transmission channel" means one or 

Let [r] denote the greatest integer which is less than or equal processing and/or radio propagation 

to real number r. Since the function q( ) is zero for negative stages between the transmitted baseband signal x(t) and the 

values of the argument and one for values of the argument received baseband signal y(t). In the preferred embodiment, 

greater than 5T, we can develop a more precise expression the term "transmission channel" includes all of the signal 

for the phase fhnction, namely ^ processing and radio propagation stages between the trans- 
mitted baseband signal x(t) and the received baseband signal 

n"^ n ^ y(Oj l^^s the transmission channel preferably includes the 

^0= ^X^* ^ 2 2j ''^^t'-^^^ RF Transmitter 130, the Radio Channel 140, and the RF 

Receiver 150. 

25 

Since, in the context of GSM, one desires to obtain a 

where n=[t/T], or in other words t=oT+>. where 0^>v<T The discrete-time linear model for the transmission channel 

current bit d and its four immediate predece^ois are stiU ^^^^^^ ^j^^ transmitted baseband 

mteracting with the filter q( ) at tune l-nT+>v, however bits , , • 5 u u j /.\ *u i • i 

which excited the filter at times kT for k^n-5 have f^lly '^^'""f ^^^^^^"^ ^hc analog signals must 

excited the filter q( ), and the re.sponse due to these bits has 30 be sampled. In the system and method of the present 

fully developed to amplitude 1. invention it is assumed that the samples of the received 

Next, in-phase and quadramre signals I(t) and Q(t) are baseband y(t) are obtained with arbitrary (but constant) 

generated by evaluating the cosine and sine respectively of phase. As mentioned above, one benefit of the present 

the phase function (j)(l); these evaluations are performed in invention is that the A/D converter 155 is allowed to 

Cosine and Sine Evaluation Units 350 and 355 respectively. 35 free-run, without external timing adjustment. Furthermore 

The in-phase and quadrature signals together define complex the samples of the transmitted baseband x(t) are obtained 

baseband signal x(t)=I(t)+jQ(t)=e''*<'\ The output 125A of during the training interval, i.e. when the training signal is 

the Cosine Evaluation Unit 350 carries the in-phase signal being asserted. 

I(t), and the output 125B of the Sine Evaluation Unit 355 _. * * ■ i 

carries the quadrature signal 0(t). Together, outputs 125A Given a discrete-time hnear system with impulse respo^^^ 

and 125B comprise the output 125 (referenced in FIG. 1) of ^^^^^^^ h(n), random input sequence X(n), and output 

the GMSK Modulator 120 sequence v(n), then the output sequence can be expressed as 

Observe that the signal x(t) is completely determined by ^ convolution sum: 
the bit stream a^ In particular, because the bits of the 
training sequence, bits a^j through agg, are system defined, ^ 
the signal x(i) is defined up to a multiplicative constant for = h ''(*)^(« " *) 
times 65T^t§87T. The last data bit immediately preceding 
the training sequence, bit a^o, stimulates the phase genera- 
tion filter q( ) at time 60T and requires duration 5T for its t ,v tu^ i« tu ; ,,1,-^ tx 

. 1^11 J 1 T-u- 1 • .L 1 L J ^r-j- where L is the length oi the impulse response n(n). lo 

output to fully develop. This explams the lower bound 65 r. i , ^ . • i r ^if ^ • ■ 

Tn- c . J * r* *L * • • u-4 * _4 calculate such an impulse response tor the transmission 

The first data bit after the trainmg sequence, bit a^?, starts . , . * i_ i r*i_ . j 

stimulating the filter q( ) and contributing to the phase at ^ necessary to choose samples of the transmitted 

time 87T Between times 65T and 87T the phase frinction t**^""!?, ^^g"^* .^W ^ ^.^^/^ "»P"^ """"^f^' 

and thus x(t) is completely determined by the training Naturally one is constrained to use the received samples 

sequence bits and is thus complete known to the receiver. the A/D converter as the output v(n) for the model. 

The transmitted baseband signal x(t), during this period 55 However since the phase used to obtain the received samples 

where the training sequence bits are stimulating the GMSK ^as not controUed, it is not known apriori how to sample the 

modulator, is referred to as the training signal. transmitted trammg signal to correspond with the received 

The in-phase and quadrature signals are fed to the RF samples. In other words, since the phase ofi&et of the 

transmitter 130. The RF Transmitter 130 images the com- received samples is not known, it is not known which 

plex baseband signal x(t) onto an RF carrier for transmission 60 samples of the transmitted training signal should be used in 

through the atmosphere, correlation with the received samples. 

The RF Receiver 150 accepts a corresponding signal from In the system and method of the present invention, a 

the atmosphere, and down -converts this RF signal to gen- search process is employed, whereby the transmitted train- 

erate a baseband signal y(t), which is a distorted version of ing signal x(t) is sampled with a plurality of phases which 

the transmitted baseband x(t). The distortions include chan- 65 vary incrementally, thereby generating a plurality of sample 

nel fading, multipath, noise, and transmitter and receiver sets (of the transmitted training signal). In other words, the 

characteristics. present invention uses a plurality of sample sets of the 
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transmiued iraining signal in performing the channel esti- 
mation. A candidate impulse response vector is then calcu- 
lated for each sample set, and the "best" candidate impulse rT"TTT"I ^ 

. , , , , ^ . for(t=0, k<P, k++) 

response vector is used as the mipulse response for the | 
transmission chamiel. 5 

A = -T 

To illustrate the process of calculating an impulse P 
response, assume that a phase A is chosen and the training foi(n-o, n<26, n-n-) 



signal samples are obtained with this phase, namely X(n)= xjn] - x(nT + 2T + a) 
x(nT+A) where n is an integer index which runs through 26 } 
consecutive values which correspond to the training ^ 



sequence bits. This involves evaluating the known math- 
ematical expression x(t) for the transmitted training signal at Since the functional expression for the training signal, x{t), 
the indicated times during the training interval. Now the is defined by the GSM specification, preferably the sample 
estimation problem specializes to finding an impulse ^5 sets are calculated and stored in the GSM receiver prior to 
response h(n) which achieves operation. 

Observe that the cross-correlation vector V, needs to be 
calculated only once for a received signal burst Furthermore 
- yhik)X{n-k). tbe correlation matrices above depend only on the known 

t=o symmetrized training sequence and the known mathematical 

^ expression for the transmitted training signal. Thus, these 
correlation matrices may be calculated pre-opcrationally, i.e. 
Next, the system and method of the present invention at compile time of software realizing the method, or at 
performs a cross-correlation between the received samples design time of hardware reaUzing the system. In the pre- 
y„ and the known symmetrized bit stream {bo,bi,b2,b3. . . . , f^ned embodiment of the present invention the inverted 
bi47}, and a cross-correlation between the training signal 25 cross-correlation matrices Mj^ are calculated and stored 
samples and the symmetrized bit stream. In view of the last pre-operationally. Thus upon receiving a signal burst, the 
equation for y„, this cross-correlation can be expressed in impulse response vectors hj^ are calculated as illustrated by 
terms of the mathematical expectations as the following code-fragment: 
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V = Corrclatcl 
for(t=0. k<P, k++) 
{ 



wherein x is an integer delay index. The cross-correlations K = M],v 

are calculated by means of time-averages over the index n ^ 

(which are explicitly specified later in this discussion). Since 

the impulse response constitutes L unknown cocflBcients, wherein the function Correlatel summarizes the calculation 

one must generate a system of L equations to determine the of vector V as described above, wherein P is the number 

impulse response. This is achieved by calculating the indi- distinct phases A, and thus the number of distinct sample sets 

cated cross-correlations for L distinct values of the delay (of the transmitted training signal) generated using the 

index x. This involves calculating L cross-correlations distinct phases. 

between the symmetrized bit stream and the received values. Furthermore, in the preferred embodiment of the present 
and cross-correlations between the symmetrized bit invention the sample sets Xo,Xj, X2, . . . ,X^_i are con- 
stream and the training signal samples (since these correla- volved with their respective impulse response vectors ho,hj, 
tions depend on the index k as well as x). In the preferred h2, . . . ,h^_j to generate a corresponding collection of output 
embodiment of this invention the linear system is solved by vectors YojY^jYi, . . . ,Yp_i as illustrated by the following 
matrix inversion. Namely, if we define the matrix R of code fragment: 
dimension LxL whose components R(x,k) are given by 
R(x,k)=E[b„X(n+x-k)], and the vector V of dimension Lxl 

whose components V(x)=E[b„y„+J, then the unknown 

impulse response may be calculated as the matrix product of for(k-o, k<P, k++) 

By repeating this procedure for a plurality of P distinct } 
phases which incrementally vary from zero to a full sjTnbol 55 

period, one generates a plurality of sample sets ^he preferred embodiment of the present invention 

Xo,Xi,X2, . . . , Xp_i firom the training signal, a plurality of sample sets from the training signal Xo,Xi,X2, . - . ,X^_, 

inverted cross-correlation matrices Mo,Mi,M2, • - - ,Mp_„ j^^st ^e stored, however they are preferably calculated 

and a plurality of impulse response vectors ho,hj,h2, . . . pre-operationally. 

y^p-v 60 Finally, each of the output vectors Yo,Yi,Y2, . . . ,Yp_j arc 

Referring now to FIG. 4, an Ulustration of the plurality of compared with a vector composed from the received values 

samples sets X^X^jX^, . . . ,Xp_i and their means of Y={yo,yj,y2,y3 • ■ • tViAi}- In the preferred embodiment of 

generation is presented, for the special case of Pt=4. In the the invention, the comparison is performed by calculating 

preferred embodiment of the present invention, the sample the Euclidean (square-norm) length of the vector differ- 
sets from the transmitted training signal Xq, Xj, X2, . . . , 65 ences. The output vector Y^^, which is closest to Y, specifies 

Xp_i are generated as illustrated in the follow pseudo-code the channel vector h^^^ which is chosen to serve as the 

fragment. impulse response for the transmission channel. This best 
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channel vector is supplied to the \^terbi Equalizer 170. This 
process is illustrated by the following pseudo-code frag- 
ment: 



MmError^LargcNumbcr 
for(k-0, k<P, k-i^) 
{ 

15 

Eiror ^ I^* (n + 4) - Y(n +• 67)1^ 

i£CErior < MinError) 
{ 

MinError « Error 
Minlndex = k 

} 



} 
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wherein the symbol |z| denotes taking the absolute value of 
a complex number. In an alternate embodiment the errors are 
calculated using the sum of absolute values instead of the 20 
sum of squares as given in the preferred embodiment. 

Referring now to FIG, 7, a block diagram of a digital 
signal processing system, in which channel estimation is 
performed according to the preferred embodiment of the 
present invention, is shown. The digital signal processing 25 
system comprises a digital signal processing (DSP) unit 700, 
a Data Store 710, a Sample Buffer 720, and an A/D converter 
155. The DSP unit 700 is coupled to the Data Store 710 and 
the Sample Buffer 720. The A/D converter 155 is coupled to 
the Sample Buffer 720. The A/D converter 155 receives a 
baseband analog signal from the RF receiver 150 (of FIG. 1), 
and samples the baseband signal at a rate of 270833 samples 
per second. As mentioned above, the A/D converter 155 is 
allowed to free-run, Le. to continuously acquire samples of 
the received baseband signal with unknown (but constant) 35 
phase. The A/D converter 155 writes the digital samples to 
the Sample Buffer 720. In the preferred embodiment of the 
invention, the Data Store 710 comprises Read Only Memory 
(ROM), Random Access Memory (RAM), FLASH memory 
or any desirable combination thereof. Tht Data Store 710 aq 
stores: executable code for the DSP unit 700; data for the 
DSP unit 7(M) to utilize in performing channel estimation; 
and the intermediate results of computations performed by 
the DSP unit 700. In the preferred embodiment of the 
invention, the Sample Buffer 720 is a Random Access 45 
Memory (RAM). 

The DSP unit 700 reads code and data from the Data Store 
710, and/or the sample buffer 720, and performs the channel 
estimation method. In so doing, the DSP unit 700 read and 
writes intermediate results of computations from and to so 
(respectively) the Data Store 710. Preferably the DSP unit 
700 is a specialized digital signal processing device. In other 
embodiments of the invention, the DSP unit 700 could be a 
general purpose processor, programmable logic, or a micro- 
controller. 55 

The Channel Estimator 160 (of FIG. 1) comprises the 
DSP unit 700, the Data Store 710, and the Sample Buffer 
720. 

The sample buffer 720 is designed to store at least 148 
successive samples of a received baseband signal y(t). Since 60 
the received baseband signal is a complex signal, i.e. com- 
prising in-phase and quadrature components, each sample of 
the received baseband signal requires two storage locations. 
Thus the sample buffer 720 comprises an array of at least 
296 storage locations. 

The DSP unit 700 monitors the power of the received 
samples in the Sample Buffer 720 to determine the begin- 
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ning of a signal burst. The power envelope of GSM signal 
bursts conform to the standard established in ETSI GSM 
05,05, European Digital Cellular Telecommunication Sys- 
tem (Phase 2), Radio Transmission and Reception, p.34, 
September 1995. The system and method of the present 
invention operates on a sequence of 148 successive received 
samples, obtained during the assertion of the signal burst. 
Since the signal burst has time duration equal to 148 bit 
intervals, it is necessary to obtain 148 successive samples of 
the signal burst, one for each bit interval. However since the 
A/D controller is free-running, without phase adjustment, 
the samples are not necessarily obtamed at times which are 
optimal with respect to the GMSK modulation scheme. 
(Recall that in the context of GSM reception, inter-symbol 
interference is minimized when samples are obtained 2.5 bit 
intervals after the bit assertion times). The Channel Estima- 
tor 160 of the present invention allows a GSM receiver to 
perform channel estimation is a way which overcomes this 
difficulty of non-optimal sample timing. Thus the A/D 
converter 155 is able to achieve a simpler and more efficient 
form, not needing the extra phase adjustment logic and 
control inputs, than in prior art receivers. 

Hereafter, these 148 successive samples, obtained during 
the signal burst, and stored in the Sample Buffer 720, will be 
denoted by the vector Y. The vector Y has length 148, and 
has complex entries. The complex entries of Y comprise 
successive samples of the complex received signal y(t). 

The Data Store 710 stores, among other things, a collec- 
tion of twelve inverted correlation matrices Mo,Mi,M2, . . , 
Mj^ (same as those described above) which are used in the 
present invention to calculate twelve candidate channel 
vectors, of which one (of these channel vectors) will be 
chosen to be the most appropriate model for the transmission 
channel. Each matrix has dimension 5x5. Furthermore the 
entries of each matrix are complex. In the preferred embodi- 
ment of the invention, the correlation matrices are computed 
and stored prior to operation in conjunction with a GSM 
receiver. The method used to compute these matrices is 
explained in detail below. 

The Data Store 710 stores, among other things, the 26 
successive bits of a GSM training sequence given by the 
values {0,0,1,0,0,1,0,1,1,1,0,0,0,0,1,0,0,0,1,0,0,1,0,1,1,1}. 
The values of the GSM training sequence are correlated with 
the received samples comprising vector Y. The result of this 
correlation is used along with the correlation matrices 
MpjMjjM^, . . . Mil calculate the twelve candidate 
channel vectors. 

The Data Store 710 stores, among other things, a collec- 
tion of twelve data vectors Xo,Xi, X2, . . . ,Xj^ which are 
used along with the candidate channel vectors to decide 
which channel vector is best, i.e. most closely approximates 
the transmission channel. Each data vector has length 26 and 
has complex entries The following pseudo-code illustrates 
the method used to generate the data vectors X^^. In the 
preferred embodiment of the invention, these data vectors 
Xjfc are computed and stored in Data Store 710 prior to 
operation of the invention is conjunction with a GSM 
receiver. 



for(k-0, k<12, k++) 
{ 

12 

for(a=0, D<26, n++) 
{ 
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Xjn] - x(nT + 2T + A) 



-continued i " 



wherein the indices x and 1 take the values zero through four. 
^ . , ^, - _ , In the preferred embodiment of ±e invention the correlation 

TTius each data vector X, comprises a sequence of 26 ^^^-^^^ calculated and stored pre-operationally. 
samples of the ftinction x(0 obtamed with a ph^^ A which ^ (comprising FIGS. 5A and 5B), 

depends on the index The funcUon x(t) is defined as the lo ^ ^^^^^^^ ^^^^^ ^^^^^^ i^^^^^^^^ ^^^^ 
baseband complex envelope which arises from a GMSK ^^^^^ ^^^^^ explained in conjunction with 

modulator, when the modulator is excited by the encoded ^^^^y^ diagram of FIG. 7. As has been mentioned, the 

bits of the GSM training sequence. Therefore x(t)=e'*^^, present invention comprises a method and system for esti- 
where the function mating a set of parameters for the transmission channel. The 

method involves receiving a signal burst from the transmis- 
" sion channel (step 500 of flG. 5A). Since the transmission 

^'^ = IZj " channel includes, among other things, the RF Receiver 150, 

receiving a signal burst from the transmission channel 
means receiving the baseband output of the RF Receiver 
The values dj, are obtained by symmetrizing and diffcren- ^ 150. [In an alternate embodiment, the transmission channel 
tially encoding the bits a^, of the GSM training sequence. The includes only the Radio Channel 140 (of FIG. 1), in which 
training bits a, are symmetrized according to the equation receiving a signal burst from the transmission channel 

b,=2a;,-l. -nic resulting symmetrized bits b, take the value "^^^"^ ^'^^^ the atmosphere ] 

1 or -1. Denote the values of the symmetrized training 25 N^^^t, the signal burst is sampled (step 510 of FIG. 5A). 
sequence by the vector B={bo,b„b3, . . . .b^,}. Tlie symme- ^Jie sampling is performed by A/D converter 155 (of FIG. 
trized bits are then differentially encoded according to the 7). llie samplmg is performed so as to obtain one sample m 
. - .u * ^ A each of a plurality of bit mtervals. Ine A/D converter 155 

equation di.-bt b«, 1, wherein the multiphcation IS a standard . .1 j- / 1 r • 1 i_ . ■ * c 1 

, . . . . , . J J . wntes the digital samples of the signal burst mto Sample 

mieger multiplication, wherein b_, is detined to be -1. ^^^^^ ^^0. The digital samples are referred to as received 

The function q(t) is defined as ^° values. Tlic received values, stored in the Sample Buffer 

720, are represented by the data structure 522 (of FIG. 5B). 
r(s)cis, ^ It is noted that, in the preferred embodiment of this 

invention, the plurality of bit intervals are contiguous (i.e. 
25 side-by-side without dapping bit intervals). However other 

1. /*\* • .1- 1 f /A J /A T-L embodiments are easily conceived where samples are 

where r(t)=g{t)*p(t),Le, the convolution of g(l) and p(t). The , . - . r ^ 1/ c i.-* - : i 

- fV. . . , 1 obtained from a plurality of non-contiguous bit mtervals. 

function p(t) is a unit ampLtude rectangular pulse with „ , . „"L ^„ . r 

turn-on interval given by OSigT. lUe function g(t) is given Referrmg now to FIG. SB, the process e^tmating 

. o J / o parameters for the transmission channel (step 520 01 rIG. 

^ 40 5 A) involves the following steps. Here, it is presumed that 

the stored sets of training signal values, vectors Xq^X^, 

- 2] X2, . . . , from the discussion above, -have been 

ln2 ' J' pre-calculated and stored in Data Store 710. These stored 

sets of training signal values are represented by data struc- 

45 ture 545 (in FIG. 5B). Furthermore, it is presumed that the 

wherein T is the symbol interval (3.69231 ^ec), and the inverted correlation matrices, matrices MoMiM2> ■ • • Mu 

time -bandwidth product BT equals 0.3. from the discussion above, have also been pre-calculated 

Recall that the Data Store 710 also stores a collection of and stored in Data Store 710. The stored matrices Mo,Mj, 

twelve correlation matrices Mo,Mi,M2. . . . ,Mii. The M^, ... ,Mii are represented by data structure 535 (of FIG. 

following code fragment describes how to generate the M^ 50 5B). Additionally, it is presumed that the values of the GSM 

juatrices training sequence {0,0,1,0,0,1,0,1,1,1,0,0,0,0,1,0,0,0,1,0,0, 

1,0,1,1,1} are stored in Data Store 710. This stored GSM 
training sequence is represented by data structure 527 (of 

FIG. 5B). 

for(k=o, k<i2. k++) 55 pi^t^ the method calculates a cross-correlation (step 525 

» lCorrelate2(X» B)]"^ ^) between the received values and the stored 

} training sequence, to produce a cross-correlation vector V, 

the same vector V as discussed above. To implement this 

calculation, the DSP unit 700 reads the vector Y of received 

where the superscript -1 denotes matrix inversion, and values from the sample buffer 720, and also reads the 

Correlate2( ,) is a function which accepts the data vector X^ training sequence from the Data Store 720, The DSP unit 

and the symmetrized training sequence vector B as input. 700 symmetrizes the training sequence according to the 

The Correlate2(,) function returns a 5x5 matrix whose equation bjt«2at-l, wherein a^^ represents the bits of the 

entries arc cross-correlations between the data vector Xj^ and training sequence. Then the DSP unit 700 performs a cross- 

the training sequence B. In particular, the (-c,l) entry, i.e. the 65 correlation between the received values and the symme - 

entry occupying the row and 1''' column, of the returned trized training sequence, and thereby generates a vector V of 

matrix is given by dimension 5x1, whose x''' component is given by 
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Next, the method calculates a plurality of channel vectors 
ho,hj,h2, . . . ,hii (step 530 of FIG. 5B), the same vectors 
ho.hithj, . . . hn 3S discussed above, using the stored 
matrices M^jM-^Mzy . ■ ■ ,Mji in data structure 535 (of FIG. 
5B). To calculate the channel vectors, the DSP unit 700 reads 
each of the matrices M;^ fr^na Data Store 710, and multiplies 
them one by one with vector V to generate a collection of 
twelve estimated channel vectors ho,hj,h2, . . . ,h^j. The 
following code fragment illustrates this calculation. 



for(k«0, k<U, k++) 
{ 

} 
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The juxtaposition of Mj^ and V is meant to imply matrix 
multiplication. 

In step 540 (of FIG. 5B), the method further involves 
producing a plurality of output vectors YojYijY^jYg, . . . , 
Yii. These output vectors are the same as the output vectors 
Yo,Yj,Y2,Y3, . . . ,Yjj explained above. One output vector 
is produced for each channel vector. To produce the output 
vectors, the DSP unit 700 reads each set of training signal 
values, vectors Xf^ from the Data Store 710 and convolves 
each with the corresponding channel vector, thereby gener- 
ating a collection of twelve estimated output vectors Yq,Y^, 
YjjYg, - . . ,Yii. The following code fragment illustrates this 
calculation. 



for(k=0, k<12, k++) 
{ 

} 



Emjii =£[Y*(n + 4)-Y(n + 67)f 
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The symbol * denotes convolution. It is to be noted that in 
an alternate embodiment, the output vector might be calcu- 
lated using a Discrete Fourier Transforms, since convolution 
is equivalent to multiplication in the frequency domain. 

Next, in step 550 (of FIG. 5B), the method determines an 
error between each of the output vectors and the received 
values of the signal burst. To determine these errors, the DSP 
unit 700 accesses the received values, i.e. vector Y, from the 
Sample Buffer 720, and calculates the Euclidean (square- 
norm) distance of each output vector Y^^ from the vector of 
received values Y. This calculation is illustrated by the 
following code fragment: 



for(k=0, k<12, 
{ 



45 



50 



60 



wherein Errorj^ denotes the error associated with the k^'' 
output vector Y;^* ^5 

In step 570 (of FIG. 5B), the method selects a best channel 
vector based on the determined errors, wherein the coeffi- 
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cients of the best channel vector specify the parameters for 
the transmission channel. The DSP unit 700 implements the 
selection by choosing the channel vector whose output 
vector Yjj, is closest to Y. In other words, the DSP unit 700 
chooses the channel vector whose output vector minimizes 
the square error. The following code fragment illustrates this 
calculation. 



M inEnor^LargcNumbcr 
for(k«0, k<12, k-M-) 
{ 

if(Errork < MinError) 



{ 



MinError «* Enor^ 
Minlndex - k 



} 



The index of the output vector which achieves the minimum 
error, Minlndex in the above code fragment, specifies the 
channel vector which is accepted as the best model for the 
transmission channel. 

Although the steps 530, 540, 550, and 570 of FIG. 5B 
have been described as separate loops, these steps can be 
equivalently implemented in one loop as illustrated in FIG. 
5C, and as in the following code fragment: 
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for(k=0, k<12, k-H-) 
{ 



Eiror = [Yt (n + 4) - Y(n + 67)]^ 



if(Enor < MinError) 
{ 

MinError = Enor 
Minlndex = k 



40 



} 



wherein the variable Error is a scalar (not a vector). In FIG. 
5C, step 530C performs the multiplication of each stored 
matrix M;^. with the cross-correlation vector V produced in 
step 525. Step 540C performs the calculation of each output 
vector Yf^ by convolving the k^^ set of training signal values 
(i.e. vector and the k''' channel vector h^^ Step 550C 
determines the square error between the k*'' output vector Y^ 
and the vector of received values Y. And Step 570C checks, 
for each new output vector Y;^) whether a minimum error has 
been achieved, and if so records the new minimum value and 
the index k. Thus when the loop has used every set of 
training signal values, the output vector which minimizes 
the square errors will have been found, and the index of this 
output vector specifies the best channel vector for the 
transmission channel. The fact that the loop is iterated so as 
to utilize all the sets of training signal values is indicated by 
the conditional branching step 560 of FIG. 5C. 

Referring now to FIG. 6, a flowchart is presented, which 
serves to facilitate the following discussion alternate 
embodiments for the channel estimation process (step 520 of 
FIG. 5A) according to the present invention. Here again, it 
is presumed that the stored sets of training signal values have 
been pre-calculated and are available. First, the method 
calculates a chacmel vector in step 630 for each of a plurality 
of sets of training signal values. The i*^ channel vector is 
calculated using the i''' set of training signal values and the 
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received values, wherein i is an integer index which runs 
through the values 0, 1, 2, ... , P-1, wherein P denotes the 
number of sets of training signal values. The i''' channel 
vector can be calculated using any system estimation 
technique, provided that the system model is supplied with 
the i*'' set of training signal values as a stimulus and the 
received values as response. 

In step 640, the parameter estimation further involves 
producing an output vector for each calculated channel 
vector, wherein the i'^ output vector is produced using the i*'' 
set of training signal values and the i channel vector. This 
involves feeding the set of training signal values to the 
system model mentioned above, wherein the parameters of 
the system model are supplied by the i'* channel vector. 

Next, in step 650, the parameter estimation method deter- 
mines an error between each of said output vectors and the 
received values of said signal burst. Any technique for 
measuring the distance between vectors may be employed in 
this step. As shown, step 660 represents the fact that steps 
630, 640, and 650 are performed for each of a pluraUty of 
sets of training signal values. 

Finally, in step 670 the parameter estimation method 
selects a best channel vector based on the determined errors, 
wherein the coefficients of the best channel vector specify 
the parameters for the transmission channel. 

Although the system and method of the present invention 
has been described in connection with the preferred embodi- 
ment and several alternate embodiments, it is not intended to 
be limited to the specific forms set forth herein, but on the 
contrary, it is intended to cover such alternatives, 
modifications, and equivalents, as can be reasonably 
included within the spirit and scope of the invention as 
defined by the appended claims. 

We claim: 

1. A method for estimating a set of parameters for a 
transmission channel in a communication system, the 
method comprising: 

receiving a signal burst from the transmission channel, 
wherein the signal burst comprises a training signal 
transmitted through the transmission channel, wherein 
the training signal encodes a sequence of binary digits 
(bits), wherein the training signal is known to a 
receiver; 

sampUng the signal burst at least once in each of a 
plurality of bit intervals and thereby producing a 
sequence of received values; 
estimating a set of parameters for the transmission 
channel, where said estimating comprises: 
calculating a channel vector for each of a plurality of 
sets of training signal values, wherein an i'* channel 
vector is calculated using an i'^ set of training signal 
values and the received values, wherein i is an 
integer index which runs through the values 0, 1, 2, 
. . . , P-1, wherein P denotes the number of sets of 
training signal values; 
producing an output vector for each calculated channel 
vector, wherein an i'* output vector is produced using 
the i^* set of training signal values and the i'* channel 
vector; 

determining an error between each of said output 
vectors and the received values of said signal burst; 

selecting a best channel vector based on said deter- 
mined errors; 

wherein coefficients of the best channel vector specify 
the parameters for the transmission channel. 

2. The method of claim 1, wherein each set of training 
signal values comprises samples of the training signal 
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obtained at substantially corresponding times in each of said 
pluraUty of bit intervals. 

3. The method of claim 2, wherein each bit interval is 
divided into P equally-spaced subintervals, further compris- 

5 ing: 

obtaining the i** set of training signal values by sampling 
the training signal at times which fall within a pre- 
defined tolerance of an i** point within an i''' subinlerval 
of said P equally-spaced subintervals; 
JO wherein i is the integer index which runs through the 
values 0, 1, 2, ... , P-1. 

4. The method of claim 3, wherein the pre-defined toler- 
ance is equal to zero. 

5. The method of claim 3, wherein the sequence of binary 
15 digits is also known to the receiver, and wherein said 

estimating the set of parameters for the transmission channel 
further comprises; 
performing a first cross-correlation between the received 
values and the sequence of binary digits to produce a 
20 correlation vector V; 

wherein said calculating the i'* channel vector comprises 
mixltiplying the correlation vector V by an i'* matrix in 
a plurality of P stored matrices. 

6. The method of claim 5, wherein the i'* matrix of said 
25 plurality of P stored matrices is calculated by performing a 

second cross-correlation between the i^'' set of training signal 
values and the sequence of binary digits, thereby producing 
a correlation matrix, and inverting the correlation matrix to 
generate an inverted matrix M^. 
30 7. The method of claim 6, further comprising: calculating 
and storing said plurality of P stored matrices prior to said 
estimating. 

8. The method of claim 1, wherein the i** output vector is 
defined as the output of a linear system suppUcd with the i''' 
35 set of training signal values as an input, wherein the behav- 
ior of the linear system is characterized by the i"* channel 
vector. 

9'. The method of claim 8, wherein the i'^ output vector is 
produced by convolving the i'^ set of training signal values 
40 and the i''' channel vector. 

10. The method of claim 1, wherein said determining the 
error comprises determining a square error between each of 
said output vectors and the received values of said signal 
burst. 

45 11. The method of claim 1, wherein said selecting the best 
channel vector based on said determined errors comprises 
selecting the channel vector whose output vector results in 
a minimum error. 

12. The method of claim 1, wherein said sampUng the 
50 signal burst comprises performing A/D conversion in which 

the sampling phase is not controlled. 

13. The method of claim 1, wherein said commimication 
system is a GSM commimication system, wherein said 
sequence of binary digits is a training sequence in a GSM 

55 traffic burst, wherein said training signal is a baseband 
complex envelop defined by the Gaussian Minimum-Shift 
Keying modulation scheme, 

14. The method of claim 1, wherein said parameters for 
the transmission channel comprise an impulse response of 

60 the transmission channel. 

15. The method of claim 1, further comprising: providing 
said parameters for the transmission channel to an equaUzer, 
wherein said equalizer uses said parameters for demodula- 
tion of the sequence of binary digits. 

65 16. The method of claim 1, further comprising: demodu- 
lating the sequence of binary digits using said parameters for 
the transmission channel. 
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17. The method of claim 1, wherein the i''' chamiel vector phirality of P stored matrices, wherein said sets of training 
comprises a candidate set of parameters for the transmission signal values and said plurality of P stored matrices are 
channel corresponding to the i*'' set of training signal values. calculated and stored in said memory prior to operational 

18. A system for estimating a set of parameters for a mode of the system. ^ 
transmission channel in a communication link, the system 5 25. The system of claim 18, wherein the i output vector 
comprising* ^ defined as the output of a linear system supplied with the 

. . , J , . • i . . r i** set of training signal values as an input, wherein the 

a receiver designed to receive a signal burs^ from the ^^^^^j^^ ^^^^ ^ characterized by the i"* 

transmission channel, wherem the signal burst com- channel vector 

prises a training signal transmitted through the trans- 26. The system of claim 25, wherein the i*^ output vector 

mission channel, wherem the trammg signal encodes a lO produced by convolving the i'^ set of training signal values 

sequence of binary digits, wherein the training signal is the i'^ channel vector. 

known to the receiver; 27. The system of claim 18, wherein said determining the 

an analog-to -digital converter configured to sample the error comprises determining a square error between each of 

signal burst at least once in each of a plurality of bit said output vectors and the received values of said signal 

intervals and thereby to produce a sequence of received burst. 

values- 28. The system of claim 18, wherein said selecting the 
a channel estimator which estimates parameters for the channel vector bi^d on said detennined errors corn- 
transmission channel; P"^f ^^<l^^& channel vector whose output vector 

, , . , results in the minimum error, 

wherein Oie channel eslimator calculates a channe ^ 29. The system of claim 18, wherein said sampUng the 

vector for each of a plurality of sets of trammg signal j ^urst is performed by an A/D converter which runs 

values, wherein an i channel vector is calculated f^^jy ^^^out controlling sampling phase, 

usmg an i set of training signal values and the 30. Jhe system of claim 18, wherein said communication 

received values, wherein i is an integer index which ij^k is a GSM communication system, wherein said 

runs through the values 0, 1, 2. ... , P-1, wherein P ^ sequence of binary digits is a training sequence in a GSM 

denotes the number of sets of training signal values; traffic burst, wherein said training signal is a baseband 

where the channel estimator produces an output vector complex envelop defined by the Gaussian Minimum-Shift 

for each calculated channel vector, wherein an i''' Keying modulation scheme. 

output vector is produced using the i''' set of training 31. The system of claim 18, wherein said parameters for 

signal values and the i'^ channel vector; the transmission channel comprise an impulse response. 

wherein the channel estimator determines an error 32. The system of claim 18, wherein said parameters for 

betweeneachof said output vectors and the received transmission channel are provided to an equalizer for 

values of said signal burst; demodulation of the sequence of binary digits, 

wherein the channel estimator selects a best channel Tbe system of claim 32, wherein said equalizer is a 

vector based on said determined errors; „ Viterbi equahzer. . . -r/. . 

wherein coefficients of the best channel vector specify The system of claim 18, wherein the 1 channel vector 

the parameters for the transmission channel. comprises a candidate set of parameters for the transmission 

19. The system of claim 18, wherein each set of training channclcorrespondmg to the 1 set of trammg signa values. 

, ^ ' r . • • - i 35. The system of claim 18, wherem said system for 

signal values comprises samples of the traimng signal ,,ti^,^^ parameters for a transmission channel is imple- 

obtained at substantiaUy corresponding times in each of said ^^^^^^ ^^^^ integrated circuits, 

plurahty of bit mtervals. 3g computer-readable storage medium comprising pro- 

20. The system of claim 19, wherein each bit mierval is gj-am instructions, wherein said program instructions arc 
divided into P equaUy-spaced subintervals; operable to implement steps for estimating a set of param- 

wherein the i''' set of training signal values is obtained by eters for a transmission channel in a communication system, 

sampling the training signal at times which fall within ^5 the steps comprising: 

a pre-defined tolerance of an i*'' point within an i*'' receiving a signal burst from the transmission charmcl, 

subinterval of said P equally-spaced subintervals; wherein the signal burst comprises a training signal 

wherein i is the integer index which runs through the transmitted rough the transmission channel, wherein 

values 0, 1 2, ... , P-1. the training signal encodes a sequence of binary digits 

21. The system of claim 20, wherein the pre-defined 50 (bits), wherein the training signal is known to a 
tolerance is equal to zero. receiver; 

22. The system of claim 20, wherein the sequence of sampling the signal burst at least once in each of a 
binary digits is also known to the receiver, and wherein said plurality of bit intervals and thereby producing a 
channel estimator is configured to perform a first cross- sequence of received values; 

correlation between the received values and the sequence of 55 estimating the set of parameters for the transmission 

binary digits to produce a correlation vector V; channel, where said estimating comprises: 

wherein said i*'' channel vector is calculated by multiply- calculating a channel vector for each of a plurality of 

ing the correlation vector V by an i"* matrix in a sets of training signal values, wherein an i^ charmel 

plurality of P stored matrices. vector is calculated using an set of training signal 

23. The system of claim 22, wherein the i''' matrix of said 60 values and the received values, wherein i is an 
plurality of P stored matrices is calculated by performing a integer index which runs through the values 0, 1, 
second cross-correlation between the i''' set of training signal 2, . . . , P-1, wherein P denotes the number of sets of 
values and the sequence of binary digits, thereby producing training signal values; 

a correlation matrix, and inverting the correlation matrix to producing an output vector for each calculated channel 

generate said i''' matrix, 65 vector, wherein an i''' output vector is produced using 

24. The system of claim 23, further comprising a memory the i''' set of training signal values and the i''' channel 
which stores said sets of training signal values and said vector; 
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determining an error between each of said output 
vectors and the received values of said signal burst; 

selecting a best channel vector based on said deter- 
mined errors; 

wherein coefficients of the best channel vector specify 5 
the parameters for the transmission channel. 

37. The mediimi of claim 36, wherein each set of training 
signal values comprises samples of the training signal 
obtained at substantially corresponding times in each of said 
plurality of bit intervals. 

38. The medium of claim 37, wherein each bit interval is 
divided into P equally-spaced subintervals, further compris- 
ing: 

obtaining the i''' set of training signal values by sampling 
the training signal at times which fall within a pre- 
defined tolerance of an i'^ point within an i'* subinterval 
of said P equally-spaced subintervals; 

wherein i is the integer index which runs through the 
values 0, 1, 2, . . . , P-1. 

39. The medium of claim 38, wherein the pre-defined 
tolerance is equal to zero. 

40. The medium of claim 38, wherein the sequence of 
binary digits is also known to the receiver, and wherein said 
estimating the set of parameters for the transmission channel 
further comprises: 

performing a first cross-correlation between the received 25 
values and the sequence of binary digits to produce a 
correlation vector V; 

wherein said calculating the i*'' channel vector comprises 
multiplying the correlation vector V by an i*'' matrix in 
a plurality of P stored matrices. 30 

41. The medium of claim 40, wherein the i''' mau-ix of 
said plurality of P stored matrices is calculated by perform- 
ing a second cross-correlation between the i^ set of training 
signal values and the sequence of binary digits, thereby 
producing a correlation matrix, and inverting the correlation 35 
matrix to generate an inverted matrix M,-. 

42. The medium of claim 41, further comprising: calcu- 
lating and storing said plurality of P stored matrices prior to 
said estimating. 

43. The medium of claim 36, wherein the i'* output vector 40 
is defined as the output of a linear system supplied with the 

i''' set of training signal values as an input, wherein the 
behavior of the linear system is characterized by the i"* 
channel vector. 

44. The medium of claim 43, wherein the i''' output vector 45 
is produced by convolving the i*'' set of training signal values 
and the i''' channel vector. 

45. The medium of claim 36, wherein said determining 
the error comprises determining a square error between each 

of said output vectors and the received values of said signal 50 
burst. 

46. The medium of claim 36, wherein said selecting the 
best chaimel vector based on said determined errors com- 
prises selecting the channel vector whose output vector 
results in a minimum error. ss 

47. The medium of claim 36, wherein said sampling the 
signal burst comprises performing A/D conversion in which 
the sampling phase is not controlled. 

48. The medium of claim 36, wherein said communica- 
tion system is a GSM communication system, wherein said 60 
sequence of binary digits is a training sequence in a GSM 
traffic burst, wherein said training signal is a baseband 
complex envelop defined by the Gatissian Minimum-Shift 
Keying modulation scheme. 

49. The medium of claim 36, wherein said parameters for 65 
the transmission channel comprise an impulse response of 
the transmission channel. 
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50. The medium of claim 36, further comprising: provid- 
ing said parameters for the transmission channel to an 
equalizer, wherein said equalizer uses said parameters for 
demodulation of the sequence of binary digits. 

51. The medium of claim 36, further comprising: 
demodulating the sequence of binary digits tising said 
parameters for the transmission channel. 

52. The mcdiimi of claim 36, wherein the i''' channel 
vector comprises a candidate set of parameters for the 
transmission channel corresponding to the i*'' set of training 
signal values. 

53. A channel estimator which estimates a set of param- 
eters for a transmission charmel in a receiving system, 
wherein the receiving system comprises: 

a receiver designed to receive a signal burst fi-om the 
transmission channel, wherein the signal burst com- 
prises a training signal transmitted through the trans- 
mission channel, wherein the training signal encodes a 
sequence of binary digits, wherein the training signal is 
known to the receiving system; 
an analog-to-digital (A/D) converter configured to sample 
the signal burst at least once in each of a plurality of bit 
intervals and thereby to produce a sequence of received 
values; 

a sample buffer configured to store said sequence of 

received values; 
a data store configured to store executable code and data, 
wherein the data comprises a plurality of P sets of 
training signal values, wherein P is an integer, wherein 
each set of ting signal values comprises samples of the 
training signal obtained at substantially corresponding 
times in each of a plurality of bit intervals; 
a digital signal processing (DSP) unit which is configured 
to read the executable code from the data store, to write 
computational results to the data store, and to read said 
received values firom the Sample Buffer; 
wherein said channel estimator comprises: 
the sample buffer; 
the data store; and 
the DSP unit; 

wherein the DSP unit estimates the parameters for the 
transmission channel by running the executable 
code, and thereby performing steps comprising: 
calculating a channel vector for each of the sets of 
training signal values, wherein an channel 
vector is calculated using an i''' set of training 
signal values and the received values, wherein i is 
an integer index which runs through the values 0, 
1, 2, . . , , P-1, wherein P denotes the number of 
sets of training signal values; 
producing an output vector for each calculated chan- 
nel vector, wherein an i'* output vector is pro- 
duced tising the i'* set of training signal values and 
the i** channel vector; 
determining an error between each of said output 
vectors and the received values of said signal 
burst; 

selecting a best channel vector based on said deter- 
mined errors; 
wherein the coefficients of the best channel vector 
specify the parameters for the transmission chan- 
nel. 
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